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Abstract: In 1981 Durhuus and Olesen (DO) showed that at infinite the eigenvalue 
density of a Wilson loop matrix W associated with a simple loop in two-dimensional Eu- 
clidean SU(A) Yang-Mills theory undergoes a phase transition at a critical size. The 
averages of det(2; - W), det{z — W)~^, and det(l -|- uW)/{l — vW) at finite N lead to 
three different smoothed out expressions, all tending to the DO singular result at infinite 
N. These smooth extensions are obtained and compared to each other. 
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1. Introduction 

So far, the large- hmit of SU(A^) YM theory has been mainly employed for qualitative 
insights. Our objective is to find a way to exploit the simplifications at large N for actual 
first-principles quantitative calculations of at least some nonperturbative quantities in the- 
ories that interact weakly at short distances and strongly at long distances. To be sure, 
significantly more work is needed in order to construct a real calculational framework; we 
are not there yet. 

In two Euclidean dimensions the eigenvalue distribution of the SU(A^) Wilson matrix 
associated with a non-selfintersecting loop undergoes a phase transition in the infinite- 
limit as the loop is dilated [1] . This phase transition has universal properties shared across 
dimensions and across analog two-dimensional models [2, 3]. Thus, a detailed understand- 
ing of the transition region in 2D is of relevance to crossovers from weakly to strongly 
interacting regimes in a wide class of models based on doubly-indexed dynamical variables 
with symmetry SU(A^). Building upon previous work [4, 5, 6], this paper presents sev- 
eral new results in this context and indicates how such results might be used to estimate 
long-distance parameters by analytical means, at least for ^ 1. 

We are focusing on the eigenvalues of the Wilson loop. The associated observables 
are three different functions /3^(^), with i = asym, sym, true, of an angular variable 9. 
At infinite A^ the three functions have identical limits: p^{9) = Poo{9)- For a specific 
critical scale, the nonnegative function poa{9) exhibits a transition at which a gap centered 
at = ±7r, present for small loops, just closes. This transition was discovered by Durhuus 
and Olesen in 1981 [1]. 

For the time being, we shall suppress the size dependence of the p%{0). p%{0) for 
i = asym, sym are obtained from the logarithmic derivative of (det'^(2; — W)) for k = 1 
and —1, respectively. One needs to take z to e*^ in a specified manner. Neither of these 
two functions p%{0) has a natural interpretation at finite A^; the interest in these functions 
mainly stems from them obeying simple partial integro/differential equations which are 
exactly integrable and already known and studied in other contexts [4, 5]. In the course of 
this paper we extend the results in [4, 5]. 

Unlike p%{0) with i = asym, sym, p%'^'^{0) literally is the eigenvalue density at finite 
A'^, (Tr (5(0 -|- ilogW)), and poses no difficulties of interpretation. It can be obtained from 
(det(l -|- uW)/{l — vW)) in the limit —u — > u — > e*^. We shall derive expressions for 
PN^^i^) t^^^ work. As anticipated in [5], we find no evidence that p%^^{6) obeys as 
simple equations as the p%{6) for £ / true were found to do. 
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This paper starts with a general description of p^^{9). We then follow with details 
for each case. First, we describe the case (. = asym, where the focus is on the loop- 
size dependence of the zeros Zj of the average characteristic polynomial. The equations 
governing these zeros were derived in [4], and here we work out the approximate solution 
for small, intermediate, and large loops. Then comes a description of the case I = sym 
and a saddle-point analysis of the integral representation found in [5]. A connection to the 
multiplicative random matrix model of [7, 8] is pointed out. We then proceed to deriving 
exact representations of p^^'^iO). As anticipated, we do not find a simple direct equation 
for p^^^^{6), but we do find a simple equation for (det(l -|- uW)/{l — vW)). From this we 
obtain a representation of p^^^^{9) by a sum. By numerically performing the sum, p^^''%e) 
can be evaluated to any desired accuracy. Further, we obtain an integral representation 
for p%^^{6) which is useful for setting up the expansion of p^™^{6). We carry out 
the saddle-point analysis which is the starting point of this expansion. We also show 
that one can define a natural extension to negative values of N, and in this extension 
P%^^{&) = P^l^{(^)- We follow this by a numerically aided study of the relations between 
the three p%{9)- We compare numerically the densities p%'^'^{6) and /0^™(6') at the same 
areas. As p^^^{0) is given by a sum of 5-functions, its comparison to another p^j^{0) is less 
direct. We conjecture, and check numerically, that the location of the N peaks in 
are close to the matching zeros 6j = —i log zj of the average characteristic polynomial. By 
"close" we mean that for large N the distance between a 6j and the matching peak vanishes 
faster than the distance between that peak and its adjacent valley. 

In order to indicate how this paper fits into a larger research plan, we finish with a 
sketch of the bigger motivating picture. 

2. Three "densities" p^{0) and how they compare 
2.1 Different definitions of dimensionless area 

The dimensionless area variable has to take a slightly different form for the average of the 
characteristic polynomial and the average of its inverse to obey equations that look simple. 
We define 

t = Ag'^N, (2.1) 

where A is the area enclosed by the Wilson loop, g is the YM coupling, and the gauge 
group is SU(A^). A = g'^N is the standard 't Hooft coupling, and t makes it dimensionless. 
This t appears in p%^°{9,t). 

The average characteristic polynomial generates the expectation values of the charac- 
ters of all antisymmetric representations of the Wilson loop matrix. The appropriate area 
variable in this case is denoted by r, with [5] 

r = t(l + ^). (2.2) 

Thus, when p^^™(0,r) is compared to p^^^^{9,t), the correction in t relative to r has 
to be taken into account. 
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The average of the inverse of the characteristic polynomial generates the expectation 
values of the characters of all symmetric representations of the Wilson loop matrix. The 
appropriate area variable in this case is denoted by T, with [5] 

T = t(l-^). (2.3) 

Thus, when /9^™(0,T) is compared to p^^^^{6,t), the correction in t relative to T has 
to be taken into account. 

2.2 Averaging over the S\J{N) Wilson loop matrix W 

The probability density for W is given by the heat kernel (see for example [9] and original 
references therein) 

VN{W,t) = J2drXr{W)e-^^^^'-^ (2.4) 

r 

with t = XA, where the sum over r is over all distinct irreducible representations of SU(A^) 
with dr denoting the dimension of r and C2(r) denoting the value of the quadratic Casimir 
on r. Xr(yV) is the character of W in the representation r and is normalized by Xr(l) = d^. 
Averages over W at fixed t are given by 

{0{W)) = j dWVN{W,t)0{W) , (2.5) 

where dW is the Haar measure on SU(A^) normalized by / dW = 1. Note that we have 
/ dWVN{W,t) = 1. Any class function can be averaged when expanded in characters using 
character orthogonality. 

Because in the sum over r in (2.4) each representation is accompanied by its complex 
conjugate representation, it is easy to see that 

{0{W)) = {0{W^)) = {0{W*)) , (2.6) 

implying identities relating {det{z — W)), (det(2; - Vl^)"^), and {det{l + uW) / {1 — vW)) to 
the same objects with z ^ 1/z, z ^ z* , u,v ^ l/u, 1/v, and u,v — > u*,v*, respectively. 

2.3 General features of the p%j{9) 

The are real on the unit circle parametrized by the angle \6\ < tt, even under 9 — > —9, 
and depend on the size of the loop. All three are positive distributions in 9, normalized by 

^P%{0) = 1 . (2.7) 

/9^^™ summarizes the averages of the characters of W in all totally antisymmetric rep- 
resentations, i.e., single-column Young diagrams. /O^™ summarizes the averages of the 
characters of W in all totally symmetric representations, i.e., single-row Young diagrams. 
/9^"^ summarizes the averages of the traces of all /c- wound Wilson loops matrices, (Tr W^). 
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As we will discuss in section 5, the latter are determined by linear combinations of the av- 
erages of the characters of W in representations which we label by (p, q) and whose Young 
diagrams have the following shape: 





1 


2 






q 


1 

2 

P 





(2. 



p^^^° determines {Ti f(W)) for any function /. However, unlike p^^™ and /O^™, it has no 
information about any average of the type (Tr/(VF) Ti g(W)), where the number of trace 
factors exceeds one. In other words, p^^^° is the single eigenvalue density and, unlike p^^™ 
and p^™, contains no information about any higher-point eigenvalue correlations. 

2.4 p%''"'{e,T) 

p^^™(^, r) is constructed from the logarithmic derivative of the average of the characteristic 
polynomial, whose zeros are Zj{T) = exp(i0j(r)) with j = 0, . . . , N — 1 and — vr < 6j < vr, 

V^(^)(z,r)E=(det(z-Vr))= 1[{z-zj{t)). (2.9) 



Define 



NipW{z,T) 



j=0 



d N 
oz 2 



(2.10) 



Setting z = e we obtain 



1 



N ^ ^ y + Oi + 2Tm 



We set y to 9 zizie with e > and real 9 and define 

9.i^)(0,r) = lim</>W(e-(^±"),r). 

Finally, we define p^^™(0,r) in analogy to p^™(0,T) in [5], 

p^y^{9,T) = -2 Re \e,T) + l]=-^ [^f )(0,r) - ^i''\9,T 

E {e + 9,{t)) = (^ - 9,{t)) , 



(2.11) 



(2.12) 



N 



(2.13) 



i=o i=o 
where denotes the 27r-periodized (5-function with normalization 



d9 62A 



I 



(2.14) 
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and the last identity in (2.13) follows from the fact that p^^^{9,t) is even in 6. The 
sum over (5-functions will reproduce exactly the averages of the traces of W in all to- 
tally antisymmetric representations at arbitrary finite A^, simply by setting W equal to 
diag(e*^''('^\ 6*^^^"^^ . . . , e*^'^-i^'^)). Thus, the entire information of p^'*'™(^,r) is contained 
in the set %(t). It is obvious that given p^^^{6,t) we can reconstruct (/>(^)(z,r) and 
ip^^\z,T). The infinite- limit of p^^™(0,r) is Poo(^,r) [2]. 

In [4] it was shown that the ^j(t) are determined by a set of first-order "equations of 
motion" in r with a specific initial condition, 

1 y-cot^^. (2.15) 

The initial condition 

^j(0) = (2.16) 

is at a singular point of the differential equations. However, once one understands that 
as r grows from zero the 0j(t) spread out, the solution becomes uniquely determined. 
Throughout the evolution, the 9j never change sign. For any r > we have 

^o(r)<^i(r)<...<^^_i(T). (2.17) 

There is a Z2 symmetry pairing them, 

ON^j^iir) = -%(t) . (2.18) 

If is odd (2.18) yields 



'JV-l 
2 



(r) = 0. (2.19) 



Thus, there are [iV/2] pairs of nonzero eigenvalues of opposite signs, implying p^^™{6, r) = 

In section 3 we shall calculate the behavior of the %(t) at small, critical, and large r. 
2.5 p%'^iO,T) 

p^™(0, T) is constructed from the logarithmic derivative of the average of the inverse 
characteristic polynomial. We reproduce here the relevant formulas from [5]. Define 

#)(z,r) = (det(z-Ty)-i), (2.20) 

where -|- is for \z\ > 1 and — for \z\ < 1. Because of the negative power, one cannot 
exclude singularities at jz] = 1 (although equation (22) of [5] shows that these singularities 
are removable so that t(^^\z,T) can be continued to \z\ = 1). One should think about 
ip^\z,T) as two distinct functions. They are simply related to each other by 

V^i^)(i/z,r) = (-z)^V'f^(z,r), \z\>i. (2.21) 

We now define 
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Figure 1: Plots of p''^''\e,T) for T = 2 (left), T ^ 5 (right), and N = 3,5,10,25,50,100,250 
together with poo{0,T). 



p^'^iB, T) is given by 

p%'^{9,T) = ilim [cpf\e-'^+%T) - cP^_^\e~'^-%T)j . (2.23) 

Unlike p^^^{9,t), p^^™{9,T) is a smooth function of for any finite N and T > 0. It 
again obeys p^^™{9,T) = p^^'^{—9,T). The function is monotonic on each of the segments 
(— TT, 0) and (0, tt) with the maximum at = and the minimum at ^ = ibTr. The infinite- 
critical point is at T = 4. For T > 4, p^™(6',T) approaches poo{0-,T) by power corrections 
in 1/iV [5]. For T < 4, poo{9, T) is zero for \9\ > 9c{T), where < 9c{T) < vr and 6*^(4) = vr. 
In this interval p^™(^, T) approaches zero by corrections that are exponentially suppressed 
in A^. /9^™(0, r) has an explicit form in terms of rapidly converging infinite sums, 

sym,. rj.. , ^ p{9,T)+p*{9,T) 

Pn i^^ T) = 1 + — , (2.24) 

PiO, T) = ^^=^^yN-,)^ ^ ^ (2.25) 



Given p^^"^{9, T) with T > we can reconstruct (l)[^\z, T) and tlj^\z, T) using the Poisson 
integral, on account of the analyticity of Tp^\z, T) for \z\ < 1. 

In [5] it was also shown that t(j^_^^\z,T) for |z| > 1 has an integral representation given 



by 



V^f)(.,r) = e^J^ rdne-^"^(.e-t-e^t)-^. (2.26) 



It was pointed out there that this formula exhibited a formal relation to (det(z — W)) 
under a sign switch of N. Similar observations have been made in the past, see [10] and 
references therein. 

Equation (2.25) can be evaluated numerically for arbitrary N to any desired precision. 
In figure 1 we show how p^^^{9,T) approaches the infinite- result poo{0) of DO [1] for 
fixed T = 2 and T = 5. In addition to these numerical results, it would be useful to compute 
analytically the asymptotic expansion of p^™(0, T) in For this it is enough to expand 
ip^\z,T) in 1/A^, which is best done by starting from (2.26). The 1/A^ expansion then 
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comes from an expansion around a single saddle point. This problem will be considered in 
section 4. The saddle points turn out to be related to the position of the boundary of the 
eigenvalue domain of the random multiplicative complex matrix ensemble studied in [7, 8] . 
In section 6 we will show plots comparing p^™(0,r) to 

2.6 p%^-{e,t) 

Finally, unlike with £ = asym, sym, t) has a natural definition. If the eigenvalues 

of W are e*"-? with j = 0, 1, . . . , — 1, we define 

pT{e,t) = -Y^{62Ae-a,{W))) = -{Tt 62.(0 + ^log{W))) . (2.27) 
j 

With the help of p^^^^{6, t) we can compute the averages of a specific subset of class functions 
F(W), namely, those that can be written as 

nW) = ^^f{a,{W)). (2.28) 

j 

The obvious formula is 

^^w) = (]^E j jom62Ao - aj{w))) = J ^£mp'r{o,t). (2.29) 

p^^'^(9,t) summarizes all the information contained in the entire collection of averages of 
the type (Tr f{W)). Viewed in this way, it is analogous to p^^^'^^O, r) and p^y"^{0, T), which 
summarize all the information contained in all averages {Xr{W)), with r denoting all totally 
antisymmetric and all totally symmetric representations, respectively. The analog of the 
functions (I)^^\z,t) and (t)~^\z,T) related to p^^™{9,t) and p^^"^{9,T), respectively, in 
the present case is the average resolvent. 

Here again the + sign goes with jz] > 1 and the — sign goes with \z\ < 1. 

Using (2.6) one easily concludes that ^^■'(z, t) determines $^^(z, t) just as in the case 
of 0i^^(z,r). Clearly, p^^^'^e^t) determines <^^^\z,t) since the latter is the average of a 
single trace. It is easy to see that the opposite is true also, namely, ^[^\z,t) determines 
p^j^^^{6,t). If we use the restrictions following from (2.6), it is enough to use just ^^_^\z,t) 
for example, 

pt™^(^, t)=2 lim Re \e'^+'<^^f\e''^+\ t)] - 1 . (2.31) 

PN^^i^i t) is smooth over the circle and similar to p^^^{9, T) in this sense, but has N peaks 
adding an oscillatory modulation to the function. In some sense p^^^°{6,t) is intermediate 
between p'^^^{6,t) and p^^"^{0,t), since it can be obtained from the expectation value of 
the ratio of values of the characteristic polynomial evaluated at two different values of its 
argument. The oscillatory behavior is in this sense a remnant of the (5- function structure 
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of p^^^{9,t). For this reason we expect the peaks of to occur at locations close 

to the matching 0j(r). This expectation will be confirmed numerically in section 6. 

Unhke for i = asym, sym, explicit formulas for t) were unavailable so far. New 

formulas that apply in this case will be derived in relative detail in section 5. We shall see 
that again a symmetry under N ^ —N exists. 

3. Motion of the zeros Zj(T) as a function of r 

In this section we only consider p^-^"^{6, r) and study the zeros zj^t) = exp(i0j(r)) of the 
average characteristic polynomial for small, large, and near the critical r. 

3.1 Oj{T) for small t 

3.1.1 Approximate "equations of motion" 

From equation (2.15) we obtain 



Rescaling 



= 4? (3-2) 



yields 

^. = EE hr^- (3-3) 

The initial condition 0j(r = 0) = indicates that one can neglect to leading order in r the 
terms with n / 0, 

- E ■ (3-4) 

In this approximation periodicity under 9j — > 6j + 2tt is lost, making the approximation 
unreliable when periodicity becomes relevant. This weak-coupling feature is a recurrent 
theme in models that have compact variables and become disordered at strong couplings. 

3.1.2 Solution of the approximate equations 

Assigning dimension 1 to r we see that t] has dimension 1/2. We define 

rjj = fijV2T , (3.5) 

making the fjj variables dimensionless and therefore independent of r. They are determined 
by the equations 

= E -4-7 • (3-6) 

The solution of these equations is well known, see, e.g., [11, App. A. 6]. The fij are the 
distinct zeros of the Hermite polynomial H]\f{x), 

^7v(%) = 0, j = 0,...,N-l. (3.7) 
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3.1.3 Relation to harmonic oscillator 



In the theory of orthogonal polynomials, the zeros of orthogonal polynomials are shown 
to be the eigenvalues of the Jacobi matrix, which is the appropriately truncated matrix of 
recurrence coefficients [12, Sees. 2.4 and 2.11]. Introduce the matrix ajy, an A^-truncated 
version of the infinite dimensional annihilation operator a normalized by 



[a, a"^] = 1 . 



(3.8) 



The truncation is to the space spanned by the harmonic oscillator states (a^")"'|0) with 
j = 0,...,N-l, 

/o \/T ••• \ 

^ ••• 



The aN satisfy 



:::::: 

••• Vin^ 

\0 • • • / 



(3.9) 



(3.10) 



where P„ = |n)(?i|. 

Using the recurrence relations of the Hermite polynomials, the Jacobi matrix is found 
to be {ajM + a^^)/V2. Thus, to leading order in r the zeros of (det(2; — W)) are the same 
as the zeros of 



det 



(3.11) 



3.1.4 Largest zeros 

Of particular interest are the largest zeros in absolute magnitude. They come in a pair of 
opposite signs. Using a known formula for large [13, Eq. (6.32.5)], we have 

1.856 



VM 

giving the largest 6j as 

0m(t) = 2^ (l 



'2N 



(2iV)i/6 



+ ... 



1.856 



+ ... 



M = N -I. 



(3.12) 



(3.13) 



(2Af)2/3 

We now are in a position to estimate when r cannot be considered to be small anymore 
and the approximation first breaks down. 

In (3.1) set j = M and choose k so that 9k = —Om- We see that by keeping only the 
re = term in the sum we neglected, for example, the following potentially large term, 

^ (3.14) 



This is the point where ignoring periodicity becomes unacceptable. At ^ 1 our small-r 
approximation breaks down for 

1.856 \ _ 

(2A^)2/3 ) 



2vV 1 



(3.15) 



- 10 - 



In conclusion, the small-r approximation holds for 

V^«| (3.16) 

a N :$> 1, but extends further if N is not too large. Since we know that at infinite N there 
is a transition at r = 4, we see that the small-r approximation cannot take us all the way 
to the critical point for iV S> 1. 

The most important conclusion is that the expansion in scale for small loops yields a 
spectrum restricted to a finite arc centered at zero angle and that the boundaries of the 
arc approach their infinite-A^ limits by a leading term of order N~'^^^. This exponent is 
a well-known property of the Gaussian ensemble of Hermitian matrices, and is connected 
to universal functions constructed out of the Airy function. The Airy function is in turn 
familiar from WKB wave functions at linear turning points. The power of 3 that appears 
in the exponent of its integral representation is related to the denominator 3 in the power 
of N we just saw. 

As the scale of the loop grows, the boundaries of the arc expand, until they meet each 
other at 9 = ±7r, at which point the small-scale expansion breaks down and the exponent 
changes. 

3.2 Ojlr) for large t 

3.2.1 The eigenvalues at r = oo 

The eigenvalues expand away from zero until they stop at r = oo, at which point they 
are equally spaced and contained in the interval (— tt, tt). Throughout the expansion they 
maintain the sum rule 

N-l 

^ej{T) = 0. (3.17) 

j=0 

This determines their asymptotic limits. 



0jir = <^) = %(j-^^)^@J, j = 0,...,N-l. (3.18) 



We now prove that the above configuration is an equilibrium point in the sense that the 
r-derivatives of the ^j(r) vanish fov 6j = Qj, j = 0, . . . , N — 1. Since 



we need to show that for each j = 0, . . . , N — 1 

1 + e'^^J^®*^) 



i 1 + e^^^i~^''^ 
2N 1 - e^(«i-«fc) ^^'^^^ 



= 0. (3.20) 



Let us denote by q the A''-roots of unity. A sum over q runs over these N complex numbers. 
We need to show that 

A = y^=0. (3.21) 
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This would then imply (3.20). The above equation already implies that the LHS of (3.20) 
is independent of j. Dividing by q the numerator and denominator of the summand and 
noticing that the restriction g 7^ 1 is identical to the restriction 1/q ^ 1 for the A^-roots of 
unity q, we get A = —A = 0. 

However, we shall soon need to evaluate other sums over q, and for these a more 
general procedure is needed. This procedure, when applied to the present trivial case, goes 
as follows. Start from 



Next, 



A 



lim 



1 + xq 
1 — xq 



1 + x 
1 — x 



(3.22) 



A = lim 
= lim 



-iV + 2 x*^ 

n>0 

-N + 



1 + x 



lim 



1 



X 



N 



1+X 

1-x 







= lim 


-A + + iV 


e-»0+ 





-iV + 2A^^a;^'^ 

A:>0 



1+X 



0. (3.23) 



This again proves (3.20). Above, we observed that will be zero if n is not a multiple 

of N, and N otherwise. We need x < 1 to perform the expansion in a geometric series, but 
at the end we can take x — > 1. Similar techniques work for all other sums over q we shall 
need. 



3.2.2 Linearization of the large-r equation 

We now expand around the infinite-r solution, to see how it is approached. From the exact 
formula for {det{z—W)) in [4], we expect the approach to be exponentially rapid, with decay 
constants given by the Casimirs of the antisymmetric representations labeled by /, where 
/ = 1,...,A — 1. This is A — 1 nonzero values, not A. The missing value corresponds to a 
uniform r-independent shift in all Oj{T), which is a symmetry of the differential equation. 
This symmetry would produce a zero mode in the linearized equation, but the mode is 
eliminated by the sum rule (3.17), which depends also on the initial condition. 
To linearize we set 

e,{T) = ej + 69,{T) (3.24) 

and expand the equation of motion to linear order in 59 j. Unlike the initial condition, 
the set provides a nondegenerate configuration around which it is straightforward to 
expand. We find 

7V-1 
fc=0 



with the matrix A given by 



A 



jk 



sin^ 



2 

E 3 eVefc for k = j 



for kj^j, 

(3.26) 
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We need the eigenvalues and eigenvectors of this matrix. Note first that Ajj does not 
depend on j and is given by 

Aji = 4 V . (3.27) 

^(l-g)(l-g-i) ^ ' 

The sum over q can be performed as before leading to 

% = • (3.28) 

Hence, the matrix A has entries Aij which only depend on (i — j) mod N. Therefore, A 



has N eigenvectors (p^^^ with components (pf\ A:, / = 0, . . . , — 1, given by 



(1^ 1 ■ 7vl(N — l) . 27rZ 7 

bf = -J^e-'^N-^ e'^^ . (3.29) 



N 

The phases have been chosen for later convenience. To evaluate the action of A on an 
eigenvector cp^''^ , we need to perform sums of the type 

^'"-^EjV^^y (-0) 

The sum over q is performed as before, and one gets 

C^'^ + —— = 21{N - I) . (3.31) 

The RHS is the eigenvalue of A corresponding to the l-th eigenvector . ^ = corresponds 
to the zero mode which does not contribute to the 69 j, so we are left with A^— 1 contributing 
modes, labeled by ? = 1, . . . , A^— 1. As expected, the eigenvalues of A come out proportional 
to the quadratic Casimirs in the /-fold antisymmetric representation, given by [14] 

C2{1) = ^1{N - I) . (3.32) 

The equations of motion (2.15) have the values of the Casimirs encoded in them. 
Thus we have found that 

69,{T) = J^Ci4^e-^^(''-'K (3.33) 
1=1 

It remains to determine the coefficients C/. Since the leading asymptotic terms at large r 
correspond to Z = 1 and I = N — 1, we only need Ci . 

3.2.3 Constraints on the coefficients 

The coefficients Ci are restricted by two quite trivial exact general properties, which imply 
for the 69kiT) that 

d9kiT) = 691{t) , 59k{T) = -59N-k-i{r) . (3.34) 
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These constraints lead to 

SOkir) = ^ pisin 

with real pi and 



N-l 



1=1 



2ttI 



(A: + 1/2) 



Pi = PN-l ■ 

Every term in the sum representing 59k{T) is invariant under I ^ N — I. 

3.2.4 Leading asymptotic behavior 

Note first that we have 

N-l 



(3.35) 



(3.36) 



(3.37) 



k=Q 



which is the term proportional to z^~^ in the expansion of (2.9) in z. 

For the leading asymptotic behavior of the 9k{j) we only need pi. We can obtain pi 
from the exact result 



1 

N 



e-2iv(^-i) 



(3.38) 



Actually, we only need this result at leading order as r — > oo. To linear order in 50k, and 
keeping only the terms with I = and / = — 1 in (3.35), we have 



N 



2i 

'n 



Pi 



N-l 

E 

k=0 



27r 



(A: + 1/2) 



Performing the trivial sum over k we get 



Pi = l- 



Hence, as r ^ oo 



or, more completely. 



<5(9fc(T) - 2 sin 



2ir 



(A: + 1/2) 



e 2JV 



vr 



^fc(^)~ ^(2A: + l-A^) + 2sin 



27r 

iV 



(A: + 1/2) 



Equivalently, we can write 



2^(^-1) 



sin(efc) . 



(3.39) 

(3.40) 
(3.41) 

(3.42) 

(3.43) 



For Qk negative the correction is positive and for Qk positive the correction is negative. 
This shows that, as expected, for increasing r each eigenvalue is distancing itself from the 
origin for all k. The correction is largest for eigenvalues in the middle of the upper and 
lower half of the circle - the eigenvalues here are the last to settle into their infinite-r 
destinations. 
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3.3 Extremal Oj^r) for r ~ 4 and large N 
3.3.1 Universal zeros 

In terms of the variable y from [4], the zeros corresponding to the angles Oj{T) mod 27r are 
given by 

g7v(^(ei(T)-7r),r) =0 (3.44) 

with 

fOO 



The universal form of qNiu^T) for large N, y 0, and r 4 is obtained by replacing the 
log(2 cosh(x/2)) above by its expansion truncated at order x^, 



2 4 
X\ X X 



log cosh - = . (3.46) 

2/ 8 192 ^ ^ 

At r = 4, we have 

/oo 
d^e-f{-y)%^iog(2cosh{./2))^ (3.47) 
-oo 

The "universal zeros" yi are defined by 

/oo 
dxe-i92(^'-^s^S'*) = 0. (3.48) 
-oo 

3.3.2 Universal numerical values 

Universal zeros have been investigated in [15]. Define 

Nx^ 4 ^xyi /o^n^ 
/in , — - — = Aifiu . (3.49) 



192 



Then 

where the /i^, A; = 1,2,... are the zeros of 

F(/x) = / due'^(^^"-"'). (3.51) 



From Table 1 of [15], we have fii = 0.8221, ^2 = 2.0227, . . . Various other results concerning 
the /ifc can be found in [15]. For the extremal positive zero at r = 4 we need to look at yl, 

yl^^^,- (3.52) 

This gives, for large N, that the zero zj^Tc) that is closest to —1 with Imzj > is 

ZM « e^(-3-7^-=^/^) . (3.53) 
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4. Asymptotic expansion of p%"^(6,T) 

The aim of this section is to construct an asymptotic expansion of p^^^{9,T) in powers of 
To this end we perform a saddle-point analysis of the integral in (2.26), from which 
p^™(0,T) can be obtained via (2.22) and (2.23). It is sufficient to study only 'ijj^\z,T) 
because il)^_^\z,T) can be obtained from (2.21). 

4.1 Saddle-point analysis 

For \z\ = 1 the integrand of (2.26) has singularities on the real-n axis. We therefore set 
z = e*^"*"*^, where e > ensures that |z| > 1 but will later be taken to zero. The integrand 
of (2.26) can be written as exp(— A^/(ii)) with 

f{u) = — +\og[ze-'-2-e'-^) . (4.1) 

We now look for saddle points of the integrand in the complex-n plane, which we label by 
u = iTU {9, T), where U {0, T) = Ur{0^ T) + iUi{9, T) is a complex-valued function of 9 and 
T. The saddle-point equation turns out to be 

-TU{e,T) U{9,T) + 1/2 _ + 

U{9,T)- 1/2 ~ ■ ^ ^ 

For e = 0, this is equation (5.49) in [8] and is related to the inviscid complex Burgers 
equation via equation (5.44) there. In the present notation, the latter equation has the 
form 

Taking the absolute value of (4.2) leads to the equation 

Vf = JJr coth{TUr + e)-U'^ (4.4) 

For e = 0, this equation has been investigated previously in [8]. However, here we keep 
e > for the time being. The singularities of the integrand of (2.26) then all have Ur < 0. 
Equation (4.4) describes one or more curves in the complex-f/ plane on which the saddle 
points have to lie (for a given value of 9, the saddles are isolated points on these curves). 
In figure 2 we show typical examples for these curves for T < 4, T = 4, and T > 4, where 
e has been chosen sufficiently close to zero. (The closed contours always enclose the points 
U = 1/2 or U = —1/2. For T > 4 and larger e, the closed contour in the left half- plane 
would be missing, but right now we are not concerned with this since we are only interested 
in the limit e — > 0"*".) Analyzing (4.2) numerically we find, for all values of T, that for a 
given value of there is always one (and only one) saddle point on the closed contour in 
the right half-plane, i.e., with Ur > 0. Note that we are showing the complex-f/ plane, in 
which the original integration contour corresponds to the imaginary axis. The integration 
contour can be smoothly deformed to go through the (single) saddle point in the right 
half-plane along a path of steepest descent. No singularities are crossed since they all have 
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Figure 2: Examples of the contours in tiie complex-t7 plane described by equation (4.4) for T = 3 
(left), T = 4 (middle), and T = 5 (right). The red dashed curves are for small e > 0, while the solid 
black curves are for e = 0. For our saddle-point analysis we keep e > 0. 

Ur < 0. There are also saddle points on the contour(s) in the left half-plane (in fact, there 
are infinitely many on the open contour), but these need not be considered. 

Once the integration contour has been deformed to go through the saddle point, we can 
safely take the limit e ^ 0"*^ . Parametrizing the contour in the vicinity of the saddle point 
hy u = u + xe*^, where x is the new integration variable corresponding to the fluctuations 
around the saddle and /3 is the angle which the path of steepest descent makes with the 
real-u axis, is given, up to exponentially small corrections in N, by 



= ^ {.e^^ + ^TUie, T))^ + logsinh " ''''' + 



(4.5) 
(4.6) 



We can now expand g{x) in x. The linear order vanishes by construction. The second 
order gives a Gaussian integral over x, resulting in 



N/2 



^l_T(l/4-C/2(e,r)) 



(4.7) 



Note that the factor e *^ cannot be pulled out of the term in square brackets because 
periodicity in 6 would be lost. 

There is a potential complication. In principle, ^"(O) and therefore the denominator 
in (4.7) could be zero, which would mean that the integral over x cannot be performed in 
Gaussian approximation. For T > 4, it is straightforward to show that g"{0) is never zero. 
For r < 4, one can use (4.2) to show that ^"(0) = only for the saddle points corresponding 
to the two angles 9 = ziz6c{T) at which poo{0, T) becomes zero (see section 2.5). This means 
that for \9\ = Oc{T) the asymptotic expansion in diverges, and that it converges ever 
more slowly as |0| — > 6^ from below. 

Note that for T < 4 and Bc{T) < 16*1 < vr the function p^y"^(6',r) is exponentially 
suppressed in N . The study of the large- asymptotic behavior in this region requires 
more work. 
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4.2 Leading-order result 

Equation (4.7) is the leading order in the 1/N expansion now show 

that it leads to p^^™{9^T) — > po^{9^T) as — > oo. We first write (4.7) in the form 

llog^f )(e^^r) = I - f{u) + 0{l/N) . (4.8) 

Note that in this order we do not need the denominator in (4.7), which corresponds to 
f"{u) (or 5"(0)). Via (2.22) and using u = iTU this leads to 

)(z, T) = i(^l- ^_l^Tu ) + ^( VA^) = -iU + 0{l/N) , (4.9) 

where in the last step we have used the saddle-point equation (4.2). Equation (2.23) then 
gives 

hm p%^i9,T)=2ReU{e,T), (4.10) 

A* — *oo 

which equals po^{9,T) of DO [1, 17] since U{9,T) satisfies (4.2) (which leads to (5.51) 
below with X = U — 1/2 and v = 1/z). 

4.3 1/N correction to p^{e,T) 

Higher-order terms in the 1/A^ expansion of vr^(e^^r) can be obtained in the standard 
way by considering higher powers of x in the expansion of g{x), resulting in integrals of 
the type J"^^ dx x^'^e~^"^^^^^^'^ with n G N. However, if we are only interested in the 1/A^ 
correction to poo{9,T) the result (4.7) is already sufficient (1/A^ corrections to this result 
would give 1/N'^ corrections to poo{9,T)). Therefore we now write 

llogV^(e*^^) = I - /(n) - ^log[T/"(n)] + 0{1/N') , (4.11) 



c^f\z,T) = -^U 1 + 4 m ^^r!:/^:^' ] + 0{1/N') (4.12) 



which leads to 

V iV[l -T(l/4-T/2)]2 
and thus to 

1 r(l/4-C/2) 



p^^-(^,r)=2Re 



U 1 + 



iv[i -r(i/4-c/2)]2 



+ 0{1/N^). (4.13) 



Note that for T < 4 and \9\ — > 9c{T) (from below) the denominator of the 1/A^ term 
approaches zero, which corresponds to the complication discussed in section 4.1. Note also 
that for r < 4 and \9\ > 9^ the saddle point U{9,T) is purely imaginary so that both 
the leading order and the 1/A^ term are zero. This confirms that the above saddle-point 
analysis is not the right tool to compute finite- effects in this region. 

In figure 3 we show examples for the 1/A^ corrections to poo{9,T) for N = 10 and 
T = 2 and 5. 
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Figure 3: Examples for the l/N corrections to Poo{0,T) for iV = 10, T = 2 (left), and T = 5 
(right). Shown are the exact result for p^™(0, T) (blue dashed curve), the infinite- TV result Poo{(^, T) 
(red dotted curve), and the asymptotic expansion of p^^'^^{9,T) up to order 0{1/N) from (4.13) 
(black solid curve). We observe that the asymptotic expansion converges rapidly for small \9\ and 
more slowly for larger 16*1. 

5. The true eigenvalue density at finite N 



We now proceed to derive exact formulas for the eigenvalue density p^"'^(0,t) 
5.1 Character expansion 

To compute (2.30) we consider the ratio of determinants 

det(l + uW) 
^ ' ' ^ det(l - vW) 

with |?;| < 1 and expand it in SU(A^) characters using [5] 



(5.1) 



det(l + uW) 



N 

E 

p=0 



det(l - vW) 



oo 



q=0 



(5.2) 



where Xpi^) iXqi^)) denotes the character of in a totally antisymmetric (symmetric) 
representation whose Young diagram consists of a single column (row) with p {q) boxes. 
The trivial representation corresponds to p = {q = 0), and for SU(A^) the antisymmetric 
representation with p = N boxes is equivalent to the trivial one because of Xjv(^) ~ 
detW = 1. This yields 



N oo 

R{u, v,W) = Y,Yl ^''^'xt{W)xUW) ■ 

p=0 9=0 



(5.3) 



The task now is to decompose the tensor product p^ (8> q^ into irreducible representations. 
In general, p^ ® q^ consists of tensors with p + q indices, where the first p are antisym- 
metrized and the last q are symmetrized. To decompose into irreducible representations we 
take one index from the first p and one from the last q and either symmetrize or antisym- 
metrize this pair. There are no more symmetrization operations we can perform. Thus, 
p^ (g) q^ decomposes into two irreducible representations, except in boundary cases when 
it is already irreducible. The boundary cases are at g = or p = or p = A^. Away from 
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the boundary cases ® decomposes into two irreducible representations identified by 
Young diagrams with the top row consisting of h boxes and a left column of v boxes and 
nothing else: 



2 



h 



(5.4) 



One either has h = q and v = p + loih = q + l and v = p. (Do not confuse the v here with 
the argument of R.) The first case corresponds to an antisymmetrized pair and the second 
to a symmetrized pair. For later convenience we shall label the "hook" diagram in (5.4) by 
{v — l,h — 1), with the understanding that -y = or /i = gives the trivial representation. 
We thus have 



{p,q-l)(B{p-l,> 



(5.5) 



and for the boundary cases 

p^0O = (p- 1,0), 0®g^ = (0,g-l) 



{N-l,q) = {0,q-l). (5.6) 



Taking into account these boundary cases and suppressing the SU(A^) matrix argument 
W, we obtain 

N-l oo Af oo 

R{u, v) = l + Y,Yl ^^^''xCP.g-i) + ^''^''X{p-i,g) ■ (5.7) 



p=0 q=l 



p=l q=0 



The case p = 0, g = is excluded from the sums. Every other boundary case appears in 
exactly one of the two sums above. Every nontrivial pair has one of the two irreducible 
representations in exactly one of the sums. Now change summation indices q ^ q + 1 in 
the first sum and p ^ p + 1 in the second to obtain 



N-l oo 

R{u, v) = l + {u + v)Y,Yl ^''^'^(P,'?) • 

This makes it explicit that R = 1 ai u = —v. 

A consequence is the character expansion of TvW^ for all k. Since 

1 



(5.8) 



R{-v + e,v) = l- — + -Tr — — + 0{e' 
V V 1 - vW 



we have 



Tr 



(p,q) 



l-vW 

where the limits on the double sum are given in (5.8). Hence, taking A: > 0, 

TrW'= Yl {-irXip,q){W) . 



(5.9) 
(5.10) 

(5.11) 



p+g — fe — 1 



Obviously, Tr 1 = iV and Tr W'' = (Tr 1^'=) 



k\* 
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5.2 Performing the average 

The average over W with weight (2.4) produces, using character orthogonaHty, 

{X(p,,){W)) = d{p, g)e- ^ (5.12) 
where C{p, q) is the value of the quadratic Casimir operator in (p, q), given by [14] 

C{p, q) = {p + q+l)(^N- ^^tA±l +q-p^ (5.13) 
and the dimension of the irreducible representation labeled by {p, q) is 

d(p,,) = d-(rtd''(,)y^^^4^^^^ (5.14) 

with 

.>,= (^), +;-!). (5.15) 

5.3 Basic combinatorial identities 

The expansions of one determinant or one inverse determinant factor (i.e., setting W = 1 
and u = ^, V = ov u = 0, V = r] in (5.8)) provide the identities 

iV-l 

S^(6 = E ^"d^iPm -P) = m + 6^-' , (5.16a) 

p=0 

°° N 
S^(r?) = Yl v'd'{q){N + q) = TT—y^ , (5.16b) 

with \ri\ < 1. These will be needed to carry out the summations over p and q later. 

5.4 Factorizing the sums over p and q for the average resolvent at zero area 

Set u = —V + e. Up to corrections of order we have 

N-l oo 

This leads to 

, > 7V-1 oo 

R{v) ^ (TV — ^ ) = - E Y.^-lfv^^'e-^'^^^^'^Uip, q) , (5.18) 

^ ' p=0 q=0 

where t = A^. Note that the sum can be extended to p = because of the factor N — p 
in d{p,q). Using (2.30) and (2.31), we obtain 

p^'^^{0, t) = 1 _ 1 lim Re[vR{v)] , v = 6"^+^^ . (5.19) 
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Note that there is no need for the hmiting procedure e ^ 0^ in (5.19) if we are using the 
double sum in (5.18) for R{v), which is well-defined for \v\ = 1. 

We now introduce an integral to get rid of the denominator in (5.14) and obtain 



R{v) 



^ dp 



TV oo 



^ ^ [{-l)PvPpPd^{p){N - p)] [v'^p'id^{q){N + q)] 



p=Q q=0 



(5.20) 



This achieves factorization of the sums over p and q at t = 0. The sum in each factor can 
be performed using (5.16). 

5.5 Integral representation at any area 

The i-dependent weight factor is the exponent of a bilinear form in p and q. By a Hubbard- 
Stratonovich transformation the dependence of the exponent on p and q can be made linear, 
and then the sums over p and q are factorized for every t and can again be done exactly 
using (5.16). 

Define the complex symmetric matrix i?jv by 



B 



N 



1 + 



N 



N N 



(5.21) 



Bn has only one eigenvalue (equal to one) and is nondiagonalizable. We have deti?iv = 1 
and 



1 



N 



' N 



1 + 



1 



\ N ^ N / 

The quadratic Casimir form can be written with the help of i?Ar: 
Cip,q)=('Py BJ'P]+N(l-l,]+N(l + l-l,]q + N(l ' 



iV2 



N iV2 



(5.22) 



2 \ 



N J 



(5.23) 



Hence 



e 2JV 



C{p,q) 



N 

— e 
t 



exp 



oo J — oo 

t 



dxdy 
27r 



exp 



Yt ^^^^ 



1 



. (5.24) 



Using (5.16) we now perform the sums over p and q, 

iV2 _ 



R{v) 



t 



X 



oo ;>oo 



dxdy 



exp 



N 
'2i 



[(1 - 1/N)x'^ + (1 + l/N)y^ - 2ixy/N] 



~oo J — oo 

"1 [1 - ^;pe-^"(*/2)(l-VAf-2/A^') 

X / dp = 

[l - ^;pe»f-(*/2)(l+l/7V-2/7V2)] 



2\-,N-l 



N+1 • 



(5.25) 
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Note that because of |f | < 1 the denominator in the last hne is never zero. The integral 
over p can be done exactly, if one wishes, resulting in 



, N _in-^) dxdy 
R{v) = — e n2) I I ^ 



f Itt 



l_ye-x-{t/2)il-l/N-2/N^) ^ 
l_y(,iy-{t/2}{l+l/N-2/N2) 



■ 4 [{l-l/N)x'^ + {l+l/N)y'^-2ixy/N] 



^ y [g-x-(t/2)(l-l/7V-2/7V2) _ giy-(t/2)(l+l/JV-2/JV2)] ' (^-2^) 

The above formula was derived for |f;| < 1; this is enough for finding p%^^{9) via (5.19). 
Using symmetries of {R{u, v, W)) one can immediately write down also results for \v\ > 1. 

5.6 Making sense of negative integer AT 

Conforming to previous observations (see [10] and references therein), we extend our result 
to negative integer N. This may be of relevance to 1/2N playing the role of the viscosity 
term in Burgers' equation [4, 5] and also to approximate equations in [16]. 
We first restate the result derived earlier, 

R{u,v,N) = {R{u,v,W)) 

N-l OO 

= 1 + ^ E E — ^«^-''e-¥^(^''''^)M^(p, N)M'{q, N) , (5.27) 

p=0 q=0 P + 9 + 



where 



C{p,q,N) = — = (p + g+ 1) ( 1 — + ) , (5.28a) 



In equations (5.28) p and q still are nonnegative integers, but N is allowed to be an integer 
of arbitrary sign (with N = excluded). 

Note that for p> N, M'^{p, N) = 0. Hence, stih keeping N > 0, we can remove one 
of the restrictions on the range of p in the sum in equation (5.27), 

OO 

R{u, v,N) = l + -uPv'^e-^^^P''^^^^M^{p, N)M'^{q, N) . (5.29) 

p,q=0 + ^ + 

Observe 

C{p,q,-N) = C{q,p,N), (5.30a) 

M^(p, -N) = {-l)P+^M^{p, N) , (5.30b) 

M^{q, -N) = {-iy+^M^{q, N) . (5.30c) 
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The entire dependence on N in (5.29) is explicit, and the function R(u, v, N) remains well- 
defined for < 0, so long as the fixed parameter is positive. With iV > this leads 
to 

R{u, v,N) = l + =^ f; (-^)^(-^^)' m^(p, -N)M^{q, -Ar)e-¥(^(9,P,-^) . (5.31) 

Interchanging the dummy summation labels p and q we get 

R{u,v,N) = R{-v,-u,-N) . (5.32) 

Writing 

R{u,v,N) = 1 + ^^n{u,v,N) (5.33) 

produces 

n{u,v,N) = n{-v,-u,-N). (5.34) 
Now set u = —V. 0,{—v,v,N) is finite for > 0. We have 

n{-v,v,N) = n{-v,v,-N). (5.35) 
Q{—v,v,N) determines p%^^{9,t) via (5.19) because of Q{—v,v,N) = NR{v), i.e., 

p%''%9, t) = 1 + A lim Re [v n(-v, v, N)] , v = 6"^+^^ . (5.36) 
At this point we realize that we have defined p^^'^^{6, t) for negative integer A^, too: 

p'l'^fie, t) = l + l. lim Re [v ni-v, v, -N)] = p'^e, t) , (5.37) 

A* e-»0+ 

where in the last step we have observed (5.35). 
5.7 Large- AT asymptotics 

If one could expand p*'^'^{9,t) in A^ around A^ = 0, only even powers of A^ would enter. 
However, all one can do is an asymptotic expansion in 1/A^, and then odd powers will 
appear. Essentially, the asymptotic expansion is not in 1/A^ but rather in 1/|A^|. For 
example, for small loops there is an arc centered at ztyr where the infinite-A^ eigenvalue 
density has a gap, and there at finite A^ one has exponential suppression of the form 
exp(— I A^|k), k > — it makes no sense to drop the absolute value on A^ in the exponent. As 
another example, consider a sub-leading term that goes like cos(A^^)/|A^|. The oscillatory 
behavior of p^"*^(0,i) comes from a contribution of this type. 

We now turn to the integral representation to take the first steps in a 1/A^ expansion 
of p^^^^{9,t). Shifting integration variables x ^ x + (t/2) [l/N + and y ^ y — 

i{t/2) {l/N - 2/N'^) in (5.25), we obtain 

(5.38) 
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Figure 4: Contours of solutions of equation (5.41) in the complex-C/ plane at t = 3 (left), i = 4 
(middle), and t = 5 (right) for p = 1 (black, solid), p = 0.9 (red, dashed), p = 0.6 (green, dot- 
dashed), and p ~ 0.3 (blue, dotted). In the figures (but not in the analysis) we have taken — 1 
for simplicity. 



Since this integral representation was derived for |f | < 1, we set v = e^^ with \6\ < vr, 
e > 0, and take the limit e — > at the end. We write (5.38) as 

^ g(Ar-l)log(l-t;pe— '/2)-(7V+l)log(l-t;pe'!'-'/2) ^ ^^^^^^ 

At large A^, the integrals over x and y decouple at leading order and can be done in- 
dependently by saddle-point approximations. Let us start with the integral over y since 
it is conceptually simpler. The y-dependent coefficient of the term in the exponent in 
equation (5.39) that is proportional to — A is 

Ky) = Y^y^ + log [i - ^pe^^"^] . (5.4o) 

Substituting y = u — it/2 = it{U — 1/2) (with u = itU in analogy to section 4) results in 
exactly the same integrand that was already considered in section 4, with the replacements 
T — > t and z ^ 1/vp (with \vp\ < 1) and with an integration over u that is now along the 
line from —oo + it/2 to -|-cxo + it/2. Since there are no singularities between this line and 
the real-u axis we can change the integration path to be along the real-n (or imaginary-C/) 
axis. Now everything goes through as in section 4. The saddle-point equation reads 

U-1/2 vp' ^ ' 

which is equivalent to (4.2). In figure 4 we show the contours in the complex-C/ plane on 
which the solutions of the saddle-point equation have to lie. (For sufficiently small p we 
now encounter the case mentioned in section 4.1 where for t > 4 the closed contour in 
the left half-plane is missing.) The relevant saddle point, which we denote by yQ{0,t,p), is 
again on the closed contour in the right half-plane. For decreasing p this contour contracts, 
but this makes no difference to our analysis. The result for the y-integral is given by an 
expression similar to (4.7). 



- 25 - 



We now turn to the integral over x. The x-dependent coefficient of the term in the 
exponent in equation (5.39) that is proportional to — is 

f{x) = - log 1 - vpe-^-*/^ = -f{ix) . (5.42) 

Substituting x = —iu — t/2 = t{U — 1/2) (with u = itU) again leads to the integral 
considered in section 4 and the saddle-point equation (5.41), except that the integration is 
now along the real-C/ axis. The positions of the saddle points of the x-integral are obtained 
by rotating the saddles of the y-integral by — 7r/2 in the complex-C/ plane, i.e., Xg = —iys- 
At a saddle point we have 

r (x,) = 1 + ^ (l + ^) = f (y,) , (5.43) 



t t \ t 

and therefore the directions of steepest descent through a saddle Us and the corresponding 
saddle Xs = —iys are identical (no rotation). By analyzing the directions along which the 
phase of the integrand is constant, we find that the integration contour can always be 
deformed to go through the (single) saddle-point in the right half-plane in the direction 
of steepest descent. Depending on the parameters p, v, and t, there is either one or no 
additional saddle point on the contour(s) in the left half-plane through which we can also 
go in the direction of steepest descent. If there is such an additional saddle point, we find 
that its contribution to the integral is always exponentially suppressed in compared to 
the saddle point in the right half-plane and can therefore be dropped from the saddle-point 
analysis. In addition, there are infinitely many more saddle points on the open contour in 
the left half-plane. However, we cannot deform the integration path to go through these 
points in the direction of steepest descent and therefore do not need to include them. An 
example for the location of the saddle points and the deformation of the integration path 
is given in figure 5. To summarize, the x- integral can be approximated by the contribution 
of the single saddle point in the right half- plane, which again leads to an expression similar 
to (4.7). 

Combining the saddle-point approximations for the integrals over x and y, we find 
that, up to exponentially small corrections in A^, the integral in equation (5.39) is given by 

Ri-) = -^e-/^ [' dp I (-^] e-- , (5.44) 



t Jq 2tt \Nf"{xo)J 

where xo = xo{6,t,p) is the dominating saddle point of the x-integral. xq is a solution of 
the saddle-point equation obtained by differentiating /(x), which can be written as 



-xo-t/2 _ ^0 

and leads to 

( 1 - 7;«fi-^o-5^ = ( 

t + Xo 



Xo -I- i 



^;pe-^o-5 ) = ( ] . (5.46) 



With (5.43) we obtain 



/-(xo)(l-.pe-o-|)^ = ^+;°(^+;°) (5.47) 
V y {t + Xo) 
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Figure 5: Example for the location of the saddle points and the deformation of the integration 
path in the complex-?/ plane for t — 5 and p — 0.95. The dashed black curves (two closed, one open) 
are the curves on which all saddle points have to lie, cf. (4.4). In this example 6 ~ 3.0. On each 
of the closed curves there is one saddle point (red dot and blue dot), and on the open curve there 
arc infinitely many saddle points, but only one of them in the region shown in the plot (green dot). 
The thin solid lines are lines of constant Re /(a;) and Re/(y). The arrows point in the direction 
of increasing Re/(a::) or decreasing Ke f{y). The dashed blue curve is the integration path for the 
y-integral along the direction of steepest descent. The solid red-blue curve is the integration path 
for the x-integral along the direction of steepest descent. 

and 

R{v) = -^e-t dp (^+;°^' , e-- . (5.48) 



Differentiating equation (5.45) with respect to p leads to 



\2 



dp pt + xo{t + xo) t + xo {t + xo) ' 



which yields 



N dxn N 

R{v) = -- dp-^ = --[xo{e,t,p = i)-xo{e,t,p = Q)] . (5.50) 

tv Jq op tv 

We know from (5.45) that xo{6,t,p = 0) = 0. If we parametrize xo{9,t,p = 1) = X{9,t)t, 
where A(^, t) has to solve 

^ = iet(A+i/2) _ 1 ' (^-^^^ 
and take the limit e — > 0"*", we end up with 

^^^^^_iVA(M^ „ = e^^ (5.52) 
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Here we need to keep in mind that we have to pick the solution of equation (5.51) which 
corresponds to the dominating saddle point xq of the x-integral for \vp\ < 1. 
Using (5.19) we obtain 

lim p^^'^e.t) = l + 2Re\{9,t) , (5.53) 

which is equal to Poo{(),t) [1, 17]. Keeping higher orders in the saddle-point approximation 
(as explained in section 4.3), we can compute the asymptotic expansion of /?^"^(0,t) in 
powers of 1/N. 

5.8 A partial differential equation for the average of the ratio of characteristic 
polynomials at different arguments 

In the expression for Q{u,v,N) that follows from (5.29) a derivative with respect to t will 
bring down the Casimir factor from the exponent. Writing 

^/ = -e^+^, i; = e^-^, fN{X,Y,t)=n{u,v,N)\t=xA (5-54) 

we can reconstruct the Casimir by derivatives with respect to X and Y. All that comes in 
is the bilinear structure of the Casimir. We obtain 



dfN ^ 1 
dt 2 



1/9 V f I d \ f d 



fN . (5.55) 



One can simplify the equation hy fj\[ ^ gj\[ = In, 



dt 2 N dYdX J 



+ ^tI^^TT^ ] 9N . (5.56) 



Rescaling X NX = Z removes all explicit dependence on N in the equation. The 
equation is linear, so we are free to rescale gi\f by any power of N we find convenient. We 
define 

Gn{Z, Y, t) ^ ^gN{Z/N, Y, t) (5.57) 

and now have 

dGN 1 / 92 a2 \ 

- ' +1^]Gn. (5.58) 



dt 2\dZ'^ dYdZ ^ 

The A^-dependence of Gat will then come in only through the initial condition at t = 0. We 
proceed to find the initial condition. Similarly to (5.16) the combinatorial factors M^'^ 
have the following generating functions: 

00 

^M^{p,N)AP = N{1 + A)^-^ , (5.59a) 

p=0 

°° N 

M'[q, N)S<^ = _ . (5.59b) 

g=0 ^ ' 
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These identities go beyond (5.16) in that they hold also for negative integer A^. Using 



1 



p + g + 1 
and the fact that at t = we have 



n{u,v,N)\t=^ = ^ \ M\p,N)MS{q,N) 



we obtain 



Observing that 



we derive 



n{u,v,N)\t=o = N' / dp 
Jo 



^1 + pu) 



dr (1 + rB)^ 



n{u,v,N)\t=o 



(1 - pv)^+^ ■ 

(l + rA)^-i 
(1 + rm^+i 



N 



u + V 



1 + u 
1-v 



N 



(5.60) 

(5.61) 

(5.62) 
(5.63) 
(5.64) 



From this we now find the initial condition associated with equation (5.58), 



GN{Z,Y,t = 0) 



-NY 



1 — eN 



+Y 



1 — eN' 



(5.65) 



The partial differential equation (5.58) and the associated initial condition (5.65) admit 
arbitrary A^, no longer restricted to integers, although for noninteger periodicity in 9 is 
lost. However, periodicity in 9 was assumed when the relation between /9^"° and was 
derived. 

One can again check whether there is a symmetry under N ^ —N . The partial 
differential equation is linear and invariant under 



Z ^ -Z , Y ^ -Y , N ^ -N . 
The initial condition switches sign under this transformation. Hence, 

GN{Z,Y,t) = -G.N{-Z,-Y,t). 



(5.66) 



(5.67) 



For noninteger N there is some subtlety in defining the cuts in the initial condition so that 
the above holds. 

By Fourier /Laplace transforms one can derive integral representations, embedding 
the initial condition at t ^ 0. To get to the density p^j^^^{9,t) via (5.19) and R{v) = 
Q{—v,v,N)/N, one needs to set u = —v, which corresponds to y = at fixed Z/N = 
— e + i9, i.e., 

pt^u<=(6),t) = 1 + ^ lim Re GN{N{-e + i9),0,t) . (5.68) 

iV e— >0+ 
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Figure 6: Plots of the densities p%'"'{e,t) (red, solid) and p'^"\9,T) (blue, dashed) for t ^ 2 (left) 
and t = 5 (right), TV = 10 (top), and iV = 50 (bottom). 



At t > the limit should be smooth, but at t = one needs to generate a (5-function 
singularity in p^^^^(6, t) at 9 = mod 2tt. We first need the y — > limit of (5.65), which is 



Gn{Z,Y = 0,t = 0) 



1 _ eZ/N 1 _ e-e+ie ■ 
Expanding the denominator in a geometric series and using (5.68) yields 



(5.69) 



?,t = 0) = l + e^^£e*'=^ + e-^^X^e- 



-ike 



k=0 



k=0 



fc=— oo 



Jke 



27r62^{e) (5.70) 



as expected. p^^^^{6,t = 0) is independent of N. 



6. Comparison of the three eigenvalue densities 
6.1 p%^^{e,t) and pT'iS^T) 

If we want to compare p^^^^{6,t) and p^™(0,T) we have to take into account the 1/N 
difference between t and T, see equation (2.3). At fixed N and t, we have to compare 
p%''%e,t) and p%'^{e,T = t{l - 1/N)). The densities p*;^"^ and p^" can be obtained 
numerically by evaluating the sums in equation (5.18) and equation (2.25), respectively. 

Figure 6 shows plots of p^^'^e, t) = pT%-^, t) and p^^^iO, T) = /)5^"(-^, T) for t = 2, 
t = 5, N = 10, and = 50 in the interval < 9 < n. As stated in section 2.5, /0^™(^,T) 
decreases monotonically in that interval. The true eigenvalue density p^^"^{9,t) has A^ 
peaks (in the complete interval [— 7r,7r]) and oscillates around the nonoscillatory function 

pr{o,T). 
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Figure 7: Plots of the density p%^'^{9,t) (oscillatory red curve) together with the positions of the 
zeros of (det(e''' - W)) (vertical blue lines) for t = 2 (left) and i = 5 (right), iV = 10 (top), and 
TV = 50 (bottom). 



6.2 p%^-{e,t) and p^^'"(0,t) 

The density p^^^{9, r) is given by a sum of (^-functions, located at the zeros of the 
average characteristic polynomial, see section 2.4. Figure 7 shows that the locations of 
these zeros are close to the positions of the N peaks of p^^^^{9, t). Here we again have to 
take into account the difference in the definitions of t and r. For fixed and t, the 
peaks of p^^^'^{9,t) have to be compared to the zeros of (det(e*^ — W)) at r = t{l + l/N). 

Computing the positions of the peaks and valleys of p*''™ and the corresponding zeros 
of the average characteristic polynomial for large N shows that the difference in position 
between a peak and its matching zero vanishes faster than the difference in position between 
that peak and the next valley. This means that 



7 

scales like 



^(peak) ^{matching zero) 



^(peak) ^(next valley) 



(6.1) 



7 oc N-f" with ^ > . (6.2) 

It turns out that the value of the exponent fi depends on t and may be different in different 
parts of the spectrum, but it is always positive (for large A^). 

In the bulk of the spectrum, the difference between peak and neighboring valley scales 
like N~'^, whereas the difference between peak and matching zero scales like N~'^ for all 
t. This results in /^buik = 1- Figure 8 shows a plot of log 7, computed for the peak closest 
to = 0, as a function of log for t = 5. The line fitted through the data points has a 
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Figure 8: Plots of log 7 for the peak closest to 6* = at t = 5 (left) and for the peak closest to 
6* = TT at r = 4 (right). Data points (red) are shown together with the fitted line (blue, dashed). 



slope of -1 + 0(10-3). (The reason for choosing 9 close to is that stable fit results can 
be obtained for lower values of N.) 

For t > 4, the infinite- limit of the eigenvalue density, PooiS,t), has no gap. In this 
case the scaling behavior does not change as one goes to higher \0\, but it is necessary to 
go to large values of N to get stable fit results for /i when \9\ is close to vr. (E.g., for t = 5 
a fit at ~ 1000 results in ~ 1.04 for the extremal peak.) 

At the transition point the situation is different. From equation (3.53) we know that 
the difference between the position of the extremal zero (the zero closest to vr) and vr scales 
Hke A^-O-'^s £qj, ^ = 4_ Between N = 1800 and N = 2800, the difference between the 
extremal zero and its critical-r approximation scales roughly like N~^''^^, the difference 
between that zero and the extremal peak position scales like N~^'^^, and the difference 
between the positions of the peak and the next valley (the valley that is closer to 9 = 0) 
scales like A^-o.83_ xhis results in //critical ~ 0.28. The plot of log 7 for that case (see 
figure 8) indicates that the value of /.icriticai niight slightly increase as one goes to even 
higher values of N (which requires more computation time). 

For t < 4 there is a gap in the spectrum. In this case, the exponent /i also has different 
values at the edge and the bulk of the spectrum, but the variation is not as large as it 
is at the critical point. E.g., for t = 3 a fit between N = 1000 and N = 1500 results in 
II « 0.64 for the extremal peak. For small \9\ we again find n = //bulk = 1- Naturally, we 
expect the exact values of the various exponents of N that enter to be rational numbers 
with denominators 3 or 4 or 12 (see section 3.1.4). 



7. The bigger picture 

For definiteness consider four-dimensional (Euclidean) pure SU(A^) Yang-Mills theory. 
Focus on the string tension and think about a lattice formulation using a single-plaquette 
action 5 = '^pS{Up), where p denotes a plaquette. For a fixed A'' one can choose s such 
that no phase transitions occur for all bare real couplings go. One can define a "string 
tension", for example, by using a Creutz ratio, 

cTcrent.iL, go) - " log ^Tr W^(L, L)) (TY VF(L + 1, L + 1)) ' ^^'^^ 
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where W{Li,L2) is the Wilson loop matrix for a rectangular loop measuring Li x L2 in 
lattice units. 

cTcreutz{L,go) Can be expanded around go ^ and go — > 00. The regimes of validity 
of these two expansions are disjoint; in between there is a crossover regime and we can 
bridge it only by numerical calculation. There are extra complications around go = 00. 
Rectangular loops of the type usually used have a "roughening" nonanalyticity in go- This 
nonanalyticity is a lattice artifact. It can be avoided by choosing loops at generic angles 
with lattice planes. Then, the definition of fxcrcutz needs to be extended. All this will 
increase the complexity of the strong-coupling expansion. At the end, only a physical 
crossover separating the ranges of the weak- and strong-coupling expansion remains. We 
have no nonnumerical calculational method to bridge it. To get the continuum string 
tension in units of the perturbative scale A we need to take the continuum limit, a correlated 
limit in which go ^ and the overall lattice scale of the loop goes to infinity. This correlated 
limit preserves the crossover. 

The idea we are pursuing is to improve the above scheme in two respects. First, since 
we wish to set up a calculation in the continuum we forget about the lattice. Instead of 
thinking about crcrcutz we consider some other observable, for definiteness the extremal 
eigenvalue 6m of a Wilson loop of size A. 

For this to make sense, we need to be able to define 6m in rcnormalizcd continuum field 
theory. We hope that this can be done by first constructing a renormalized polynomial 
in z corresponding to {det{z — W)) and taking the roots of it to define 9m- While we 
have some idea how a calculation for small loops might proceed, for large loops we need 
something beyond ordinary field theory. Here we assume that an effective string model will 
describe (dct(z — W)). This model will have a dimensional parameter, the string tension, 
and will be a good description for very large loops, with corrections parametrized by more 
parameters becoming more and more important as the loop shrinks. 

To relate the string tension to A, the dimensional parameter entering the perturbation 
theory for small loops, one needs to join the two regimes over the crossover. Here is the 
point that the simplification of large N enters: At infinite the crossover for zm collapses 
into a point, and we have a phase transition. We postulate that we know that the transition 
is universal and that we know it is in the same universality class as the DO transition. 

Therefore, for N 1, the dependence of zm on intermediate scales, i.e., scales in the 
vicinity of the critical scale, is known up to a few constants. This is the ingredient that was 
missing in the lattice scenario described above. It is now possible to imagine calculating to 
some order at short, intermediate, and long scales, and sew together the three scale ranges. 
Requiring smooth matches could produce a number for the string tension in units of the 
perturbative scale A. 

There are many variations possible, zm is only one possible example of a potentially 
useful variable, zm depends on the dilation of a fixed-shaped Wilson loop, measured by a 
dimensionless variable A. As a function of A, zm will trace out a trajectory from 6 = at 
A = (one could replace this by a < Aq ^ 1) to = 7r(l — 1/2A) at A = 00. For small A, 
the perturbative scale A enters the calculations, and for large A the string tension enters. 
The two regimes are joined by the crossover. To parametrize the crossover one has to work 
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out the details of the two-dimensional case. 

In two-dimensional YM renormalization is trivial, perturbation theory is well-defined, 
and there even exists an exact string description [9]. We need to gain control over the 
crossover at large N, and then we can try to build a prototype of the calculation we 
envisage. We also need to learn enough to open the possibility of finding other interesting 
observables than 9m- This is where the present paper fits in. 
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